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Abstract. In this paper we propose a thermodynamically consistent model for superfluid-normal 
phase transition in liquid helium, accounting for variations of temperature and density. The phase 
transition is described by means of an order parameter, according to the Ginzburg-Landau theory, 
emphasizing the analogies between superfluidity and superconductivity. The normal component of 
the velocity is assumed to be compressible and the usual phase diagram of liquid helium is recov- 
ered. Moreover, the continuity equation leads to a dependence between density and temperature in 
agreement with the experimental data. 
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Introduction 

The phenomenon of superfluidity occurs mainly in liquid helium below a characteristic temperature 
0\. Above 9\, helium behaves like a conventional fluid with small viscosity. However, when the 
temperature is lowered below 8\, liquid helium undergoes a phase transition characterized by the 
ability of the liquid to flow across narrow channels without apparent friction. Helium has two stable 
isotopes 4 He and 3 He that become superfluid at low temperatures. The most common isotope is 4 He 
whose transition temperature, called the A— point, is about 2.17K. The normal phase of 4 He is called 
the He I-phase and the superfluid state is said He II. 

In this paper we propose a phenomenological model to describe the phase transition in 4 He. The 
first model to study the behavior of 4 He was the two- fluid model, suggested by Tisza ([17]) and 
developed by Landau ([12]). According this theory, when the temperature is under 6\, each particle 
of the fluid is endowed with two different excitations at the same instant: one of these is the superfluid 
velocity, denoted by v s , the other one is the normal velocity v n . The density of the fluid is the sum 
of a normal and a superfluid component 

P = Pn + Ps 

and the total current density is given by 

j = Priori + PsV s . 

If the temperature overcomes the A— point, the density p s vanishes, so that liquid helium becomes a 
normal fluid. 

The two-fluid model has been widely adopted to describe some phenomenological aspects of su- 
perfluidity, when the involved velocity of the fluid is quite small (see [16] and references therein). 
Besides the two-fluid model, some authors follow the one-fluid theory of liquid helium, which is based 
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on the extended irreversible thermodynamics (see [13, 15] for instance). They analyze the behavior 
of liquid helium II considered as a unique substance obeying a suitable Navier-Stokes equation. More 
recently in [6] the author studies the phase transition between helium I and helium II in the frame- 
work of the Ginzburg-Landau theory, by considering this passage as a second order phase transition 
and introducing a scalar variable cp as order parameter such that ip 2 represents the concentration 
of the superfluid phase. This point of view emphasizes also the analogies between superfluidity and 
superconductivity ([14, 16]). Indeed, as in the two-fluid model, the velocity of the fluid is due to a 
normal and a superfluid excitation; however the superfluid component v s is supposed to satisfy an 
evolution equation similar to the differential equation governing the motion of the superconducting 
electrons inside a superconductor ([5]). 

In our paper we consider a generalization of this model by keeping into account variations of the 
mass density of the fluid and variations of the temperature. The main assumption, distinguishing our 
model by the one proposed in [5], is that the normal component is a compressible fluid. Accordingly, 
the pressure becomes a new variable of the problem and the divergence of the normal component v n 
satisfies a constitutive equation depending on the phase variable tp. When the fluid is in the normal 
state, the evolution equation for the normal component v n reduces to the Navier-Stokes equation. 



V 


Solid 








A— line 




helium II 


i helium I 





e 



Figure 1. Phase diagram of 4 He. 



The occurrence of the pressure in the phase equation allows us to recover in the phase diagram the 
line (A— line) separating the normal from the superfluid phase (see Fig. 1). 

We prove that our model is consistent with thermodynamic principles, deducing the differential 
equation for the temperature from the energy balance law and proving that Clausius-Duhem inequality 
is satisfied. In particular we assume that the heat flux is the sum of two contributions: the first term 
is proportional to the gradient of the temperature (Fourier law) , the other one is due to the superfluid 
transition since it involves the superfluid component v s and the phase variables ip. This point of view 
is not dissimilar from models based on the extended irreversible thermodynamics, which interpret the 
superfluid velocity kind of heat flux ([13, 15]). 

In the last section of the paper we show that the differential system governing the evolution of the 
fluid can be written by means of a different set of variables, similar to the unknown fields used in the 
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context of superconductivity ([7, 18]). Such a formulation could be useful in proving some analytical 
results concerning the well-posedness of the system. 

1. The Ginzburg-Landau equation for the superfluid concentration 

As known, the passage from the normal phase to the superfluid state is a second-order phase transi- 
tion, since no latent heat is involved [14]. Therefore we propose a model to describe the phenomenon 
in the context of the Ginzburg-Landau theory. The first step is the identification of a suitable order 
parameter characterizing the state of the material. Here we introduce a scalar variable <p (phase- field), 
such that ip 2 represents the concentration of the superfluid phase. Thus the values of ip 2 are bounded 
in the interval [0, 1] with p 2 = in the normal phase and ip 2 / in the superfluid regime. The vari- 
able <p provides a measure of the internal order structure of the material, since the superfluid phase 
is considered a more "ordered" state than the normal one [12]. Accordingly, the differential equation 
governing the evolution of <p can be interpreted as a balance law on the internal order structure. For 
a general treatment of balance laws in continuous bodies with microstructure see [3, 4]. The inter- 
pretation of the Ginzburg-Landau equation as a balance equation has been proposed by Fried and 
Gurtin who introduce the notion of microforces ([10, 11]). In this context every change of the order 
parameter is related to the existence of microforces which expend power on the atomic configurations 
inside the material. A similar approach has been proposed in [5] where a balance of the internal order 
structure is postulated. The common idea of the two approaches is that, during the transition, in 
any sub-region A of the body, the power expended on the atoms by the lattice is balanced by the 
power expended across the boundary dA by the configurations external but neighboring to A and 
by the power expended by sources external to the body. We briefly recall the interpretation of the 
Ginzburg-Landau equation as a balance law, adopting the terminology of [5]. 

Let us consider a superfluid occupying a bounded subset Q C M 3 with regular boundary dQ, whose 
outward normal is denoted by n. For any sub-body A C fl, we denote by S l (A) the rate of absorption 
of the order structure per unit time, defined as 

(1.1) 5*(A) = / pkdv, 

J A 

where p is the mass density and k is the internal specific structure order. Similarly, the external order 
structure S e (A) is written in the form 

(1.2) S e {A) = [ p-nds+ / padv, 

JdA J A 

where the vector p denotes the order structure flux coming from the boundary and a is the structure 
order supply. 

Hence, the order structure balance is expressed by the equality 

(1.3) / pkdv = / p-nds+ / padv Vicfi. 

J A JdA J A 

In local form, the integral equality (1.3) leads to the equation 

(1.4) pk = V-p + pa. 

Hence, the quantity p can be assimilated to a vector stress and k and a to internal and external 
microforces distributed in the domain Q. As usual in phase transition problems, we assume a = 0. 
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where the superposed dot stands for the time derivative, m > is a suitable coefficient depending 
on the variables that induce the transition, the potentials F and G characterize the order and the 
feature of the transition and k, t are positive constants. Accordingly, the evolution equation for the 
order parameter reads 



In steady and homogeneous conditions (i.e. p = and Vp = 0) the solutions of (1.7) are the 
stationary points of the function 



Such expressions are the same used in the classical Ginzburg-Landau theory of superconductivity 
([9]). The functions F and G satisfy the following properties (see Fig. 2): 

(i) the function W admits its minimum value at <p = when m > 9\ and at p = ±p>o(m), 
< po(m) < 1, when < m < 8\; 

(ii) <po(m) — > as m — > 9\ and (fo{m) — > 1 as m — > 0. 



Figure 2. Plot of the potential W for different values of m: in the left m > 0\, in 
the middle m = 9\, in the right m < 9\. 

As a consequence the transition occurs when m = 9\ and m> 9\ identifies the normal phase. 

Notice that there is no physical distinction between the positive and negative values of (p since the 
physical quantity is ip 2 . As a consequence, we can consider p G [0, 1] as well as <p £ [—1, 1]. In the 
latter case, we have to require that the potentials F and G are even functions. 

Finally, we associate to equation (1.7) the usual Neumann boundary condition 



(1.7) 



rpp = • ( P Vp) - P 9 X F'(p) - P mG'{p). 



W(p) = 9 x F{p) + mG(p). 
A typical choice adopted for second-order phase transitions is 



(1.8) F{p) = \-\, G{fp) = 




(1.9) 



p • n\ dn = 0. 
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2. Evolution equation for the velocity 

According to Landau's viewpoint, we assume that each particle of the superfiuid exhibits a normal 
and a superfiuid excitation. Thus the velocity v is given by the superposition of such flows. In 
particular, we let 

(2.1) v = {l-p 2 )v n + p 2 v s . 

where v n , v s will be called normal and the superfiuid components, respectively. This does not mean 
that superfiuid is considered as a mixture of two fluids. Indeed, a particle of the material is endowed 
with two simultaneous excitations (normal and superfiuid) when p ^ 0, while its velocity coincides 
with v n when p = 0. However, the expression (2.1) of the velocity is not in contrast with the 
traditional theories of superfluids which assume 

pv = p n v n + p s v s , 

when p n and p s are identified respectively with p n = p(l — p 2 ) and p s = pp> 2 . 
The variables p and v are related by the continuity equation 



(2.2) /j + pV-v = 0. 

By paralleling [6], the differential equations governing the evolution of the component v s are 

(2.3) V s = -V(P S - V x v„ - pp 2 v s + V9 

(2.4) V-[pp 2 v s ] = -TK 2 pp 2 cP s 



where 4> s is a suitable scalar function referable to a "pressure" due to the superfiuid component. It 
is worth noting that (2.3) and (2.4) are similar to the evolution equations governing the motion of 
superconducting electrons [9] , emphasizing the evident analogies between superfluidity and supercon- 
ductivity. 

We associate to these equations the boundary condition 

(2.5) v s -n|ac = 0. 

A further boundary condition on v s has to be prescribed. Our model allows us to choose such a 
condition in an arbitrary way, since, as we will see in Sect. 4, the condition (2.5) is sufficient to 
ensure the vanishing of the power flux at the boundary of the domain. In the superconductivity 
model it is assumed 

(2.6) (V x v s ) x n\ dn = u, 

where a; is a known function. Therefore we can assume (2.6) by analogy with that model. 
For the velocity v n we propose the following equation 

(2.7) p(l - p 2 )ir n = V x v s + z>V(V • v„) - Vp - \ V • (Vp Vp) - pp 2 V9 + pg, 

where v is the viscosity coefficient, p is the pressure and g denotes the external force density. By 
applying the curl operator to equation (2.3), we obtain 

V x v s = —V x V x v„ - V x [pp 2 v s ]. 

A substitution into equation (2.7) leads to the generalizing Navier-Stokes equation 

p(l - p 2 )v n + V x V x v„ - z^V(V • v„) = -Vp - \ V • (Vp ® Vp) - pp 2 V6 - V x [pp 2 v s ] + pg. 
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We append to (2.7) the usual boundary condition: 
(2.8) v n | an = 0. 

3. Phase diagram 

The differential equation of the phase variable is completed by the following constitutive choice: 

m = 9 + \p + v 2 s -vl, 

Hence, (1.7) reads 

(3.1) rp0 = V • (pVp) - p0\ <p(<p 2 - 1) - p{9 + Xp + v 2 - v 2 n )<p. 

As pointed out in Sect.l, in homogeneous and steady conditions the material is in the normal state 
when m > 6\, that is when 

(3.2) 9 + Xp + v 2 s - v 2 n > 9 X 

and the transition occurs when 9 + Xp + — = 9\. Therefore the model is naturally able to 
account for the existence of a critical velocity v s (depending on the temperature and the pressure), 
above which superfluid properties disappear (see [16]). 

The regions where 9+Xp+v 2 — is respectively grater and smaller than 9\ are the stability regions 
of the normal and superfluid phase and the curve represented by the equation 9 + Xp + — = 9\ 
separates such regions. In particular, if we consider the equilibrium states, i.e. v s = v n = 0, then 
the curve is represented by the equation 

+ Xp = 9 X , 

which is a line with negative slope — 1/A. This is good approximation of the A— line shown in the 
phase diagram of liquid helium and represented in Fig.l. 

In a previous model for superfluidity ([6]) the coefficient A was supposed to vanish. This corresponds 
to approximate the A— line with a vertical line. Moreover in such a model the normal component was 
assumed to be incompressible, so that 

V • v n = 0. 

In this paper, we consider a different case, by assuming that the normal component obeys the con- 
straint 

(3.3) V • v n = A/w 
which generalizes incompressibility. 



4. Heat equation and thermodynamics 

In order to obtain the kinetic equation for the temperature, let us consider the first law of thermo- 
dynamics in the form ( [8] ) 

(4.i) P E = r^ + K s + P h, 

where E is the total energy, V % ip ^V 1 Wa are respectively the internal powers due to the order parameter 
and to the velocity v s and h stands for the rate at which the heat is absorbed by the material. 
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Hereafter we will consider some approximations of our model, valid in a neighborhood of the transition 
temperature. In particular, in the expression of the time derivative 

ip = d t <p + v • Vtp = d t p + v n • Vip + p 2 (v s - v n ) • Vp 

we neglect last term, so that p assumes the form 

<p = d t <p + v„ ■ Vp. 

As a consequence 

(4.2) Vp = (VipY + (Vv n ) T V<^, 

where the superscript T denotes the transpose of a tensor. 

Multiplying equation (3.1) by p and accounting for (3.3) and (4.2), we obtain the power balance 
related to ip, that is 

<pi _ -pe 

where the internal and external powers are given by 

d_ 
dt 



K = P: 



+ p(0 + v 2 - v 2 )pp 



(4.3) -v n -Vp + T P <p 2 + \vv n - (V<p®V<p), 

(4.4) V% = V ■ (jLp<pV<p-pv. 

Similarly, by multiplying equation (2.3) by v s + V(p s — V9, we obtain 

|v s + V<p s - V9\ 2 + v n • V x ir s + pp 2 v s ■ (v s + V<p s - V0) 
= -V • (v„ x v s + v n x V(j) s - v n x V9). 

We substitute the term V x ir s with equation (2.7) and we take (2.4) into account. Thus, we obtain 

T>i — -pe 

where the the internal and external powers due to the velocity v s are 

(4.5) V l Vs = \v s + VcP s -V0\ 2 + \pj t \v n \ 2 ~ pp 2 v n -v n + v\V -v„| 2 + v n -Vp 

— \ Vv„ • (Vp <8> Vp) + pp 2 v s ■ v s + TK 2 pp 2 (p 2 - pp 2 (v s - v n ) • V9 

(4.6) V% s = - V • [v n x v s + v„ x V^ s - v n x - w„ V ■ v n + ^ (V^ 8 V^)v n 



+pp 2 v s (j) s ] + pg ■ v n . 



It is worth noting that, as a consequence of the boundary conditions, after an integration on the 
domain Q the only contribution given by the external powers is due to the external source g. 
We assume that the total energy E is written as 

(4.7) E = ^|Vvf + 9 X - £\ + e (6) + \p\ 2 s + 1(1 - ^)v 2 , 
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where eo(0) is a function depending only on the temperature. We identify the first three terms of 

(4.7) with the internal energy 

+*<»>. 

Last two terms of (4.7) involving the normal and the superfluid components of the velocity define the 
kinetic energy T. Hence 

E = e + T. 

Substitution of (4.3), (4.5) and (4.7) into (4.1) yields 

ph = pe' (9)8- p9<pip-Tp<p 2 -\v s + X7(/) s -V6\ 2 

(4.8) -z/|V • v n | 2 - TK 2 P (p 2 (f> 2 s + W 2 (v s - v„) • V0. 

The heat equation is given by 

(4.9) ph = -V -q + pr, 

where q is the heat flux and r is the heat supply. In this framework, the heat flux is assumed to 
satisfy the constitutive equation 

(4.10) <l = -ko(0)V0-p<f?6(v 8 -v n ), 

where &o(0) denotes the thermal conductivity. Notice that, when the fluid is in the normal phase, i.e. 
ip = 0, equation (4.10) reduces to the usual Fourier law. On the contrary, in the superfluid state, the 
superfluid component of the velocity v s is related to the heat flux inside the material. 

By comparing (4.8) with (4.9), we obtain the evolution equation for the temperature, i.e. 

pe' {6)6 = p9cp(p + T P <p 2 + \ir s + V(f) s -Ve\ 2 + v\V -\r n \ 2 

(4.11) +TK 2 p^ 2 <j) 2 + V • [*o(0)V0] + V • [ W 2 (v s - v n )]0 + pr. 

Equation (4.11) is completed by the boundary condition 

q-n| an = 0. 

In view of the constitutive equation (4.10) and the boundary conditions (2.5) and (2.8), we conclude 
that 

(4.12) V0-n|an = O. 

Now we prove that our model is consistent with the second law of thermodynamics. We write the 
Clausius-Duhem inequality in the form 



rt*-v-(?) + s, 



where r\ is the entropy. We introduce the Helmholtz free energy density ^ = e — r/6 which is supposed 
to depend on (0, ip, V<p). In view of (4.1) and (4.9), we deduce 



p* ■ p n o ■ pi <: v. ■ p; ( ' q - V0. 
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Accordingly, we have 



p(d e ty + rf)6 + p[d v V - 9\ p(p 2 - 1) - 9<p]<p + p 



-rpp 2 - \v s + - V9\ 2 - v\V ■ v n | 2 - TK 2 pp 2 <\> 2 s - ^p-\V9\ 2 < 0. 

V 

In virtue of the arbitrariness of (9,p,'Vp), the functions q, ty,r/ are compatible with the second law 
of thermodynamics if v, r, ko assume non-negative values, and the free energy ^ satisfies the following 
conditions: 

d e V = -r), d v V = 9\ p{p 2 - 1) + 9p, ckj^ = ^ 



Hence, 

(4-13) * = 9 X - ^ + l -9p 2 + ^ |V^| 2 + MO), 

where ^0 depends only on the temperature. From the relation ^ = e + 9dg^ it follows that 

%(9) = tf o (0) ~ eo(0). 

Finally, we conclude this section by proving that the passage from the normal phase to the superfluid 
one is a second-order transition since no latent heat is involved. Indeed, from (4.13), it follows that 
the entropy assumes the form 

■n = -d e * = -\<p 2 -*' Q (0) 

and the latent heat L is given by 

L = 9[ v (9 x ,p (9 x ))- V (9 x ,0)], 

where po is the minimum of the function W(<p) = 9\F(p) + G(ip) and p = po, p = characterize the 
pure phases (see Fig.2). Since po = when 9 = 9\, we have that L = 0. 

5. The differential system and gauge invariance 
Collecting the equations of motion we write the system of equations: 

(5.1) rpp = \v -(pVp)- P 9 x p{p 2 -1)- p(9 + \p + v 2 s -v 2 n )p. 

(5.2) V • [pp 2 v s ] = -TK 2 pp 2 4> s 

(5.3) v s = - V0 S - V x v n - pp 2 v s + V0 

p(l - p 2 )ir n = - V x V x v„ + vV{V ■ v„) - Vp - -3 V • (Vp <g> Vp) 

(5.4) -pp 2 V9 - V x [pp 2 v & ] + pg. 

(5.5) V • v n = Xppp 

(5.6) p = -pV ■ [(1 - p 2 )v n + p 2 v s ] 

pe' (9)9 = p9(p(p + Tpp 2 + \v s + V<p s -V9\ 2 + v\V • v n | 2 

(5.7) +TK 2 pp 2 (j ) 2 s + V • [k (9)V9] + V • [pp 2 (v s - w n )]9 + pr, 
in the unknowns (p, v s , S , v n , p,p, 9. 
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As we have pointed in Section 2, the model we propose is similar to the Ginzburg-Landau model 
of superconductivity. In order to stress this analogy, following [2], we introduce the transformation: 

(ip,v s ,(f) s ) < — > 0,A,c/>) 

where 

<p = ipe~ iX , v s = A- 1 Vx, (p s = (t> + -x, 

K K 

X is an arbitrary scalar function and i denotes the imaginary unit. 

Our aim is to write equations (5.1)-(5.7) by means of the variables (ip, A, (f). Multiplying equation 
(5.1) by e tx , we obtain 

rpi) = (iT P fx+ \p^ + ^Vp-V^) e ix - MMIVf -l)+0 + - v 2 ] 

V K tx J 



(5.8) _^(|A| 2 --A-Vx + ^|V X | 2 



Dividing (5.2) by nip and substituting the expressions of v s and (p s , we deduce 

11 2 2 

rp<PX = • Vp + -^Vp ■ Vx pA ■ V<p + -i,p\7ip ■ Vx 

K K K K 

(5.9) pwV • A H — ^pipAx — TKpuxj). 

K K 

We substitute equation (5.9) into (5.8). Moreover, the identities 

W = (\7(p + i(p\7x)e ix , 

AV = (Atp + 2iV(p- V X - tp\V X \ 2 + i<P&x)e ix , 

lead to 

1 1 i _ 2i 

rpip = -^pAtp + -^Vi/j ■ \7p tpA • \7p pA • Vip 

K z K z K K 

--pipV ■ A - irnpipcf) - p^[6 x (\^\ 2 - 1) + 9 + \p - v 2 ] - pV|A| 2 . 

In addition, it is easy to prove the relations 

ir s + V(f> s = A + \74> 

<p\ s = -I^A+^-CW-V'VV') 
2k 

if 2 + kVVs = \lp\ 2 + K 2 (p 2 (j) 2 + 2K(j)ip 2 X 



\ip\ 2 + K 2 \ip\ 2 4> 2 - iK<j)(ipip* - ifnp* 

1 



\7ip®\7(p = — 2(ip*Vip + i/jVi/)*) ® (Y>*V^ + ^V^* 
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where tp* denotes the conjugate of ip. Accordingly, the equations (5.1)-(5.7) transform into 
rpip - 



1 i 1 _ 2i 

-^pAip ipA ■ Vp + -^V^j ■ Vp pA ■ Vip 



K 



— pV>v • a - impact) - pV^aGV'I 2 - !) + 6 + x p - v n] - P^\M 

K 



-V^-Vxv„| 



2k 



(y>vy>* - W>) + v# 



p(l-|Vf)v„ = -V x V x v n + zA7(V- v n ) - Vp 



1 

4k 2 



V 



1 



(V>*vy> + y>vy>*) ® (Y>*vv> + ^W) 



V-v n 
P 

pe' (o)e 



A 



-p|^| 2 V0- Vx 

p{Tpip* + 



/9|V»| 2 a + — p^vV'* - ^*W>) 



-pv • [(i - |^| 2 K - |^| 2 a + ^(VW - V*VV)] 

= -p9(lp7p* + 1p*1p) + T/9|V>| 2 + TpK 2 \lp\ 2 (f) 2 - iTKp4>(lpl(j* - Iplp*) 



+| A + V(j) - V6\ 2 + u\V ■ v n | 2 + V • [ko(9)V0] 



+v 



P \M\ n - p \^\ 2 a + ^p(^vv* - r^) 



9 + pr. 



It is worth noting that the previous equations are independent of \-> which can be chosen arbitrarily. 
This formulation allows a more direct comparison with the Ginzburg-Landau model of superconduc- 
tivity, where the choice of the variables (ip, A, cp, 9) is crucial in order to prove well-posedness results 
for the differential system. 
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